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5— ( , We study the gauge sector of Minimal Walking Technicolor, which is an SU(2) gauge theory 

with Uf — 2 flavors of Wilson fermions in the adjoint representation. Numerical simulations are 
■^L ' performed on lattices Nt x A^^, with Ns ranging from 8 to 16 and Nt — 2Ns, at fixed /3 = 2.25, 

and varying the fermion bare mass mo, so that our numerical results cover the full range of fermion 
masses from the quenched region to the chiral limit. We present results for the string tension and 
the glueball spectrum. A comparison of mesonic and gluonic observables leads to the conclusion that 
the infrared dynamics is given by an SU(2) pure Yang-Mills theory with a typical energy scale for 
the spectrum sliding to zero with the fermion mass. The typical mesonic mass scale is proportional 
to, and much larger than this gluonic scale. Our findings are compatible with a scenario in which 
the massless theory is conformal in the infrared. An analysis of the scaling of the string tension with 
the fermion mass towards the massless limit allows us to extract the chiral condensate anomalous 
dimension 7,, which is found to be 7* = 0.22 ± 0.06. 
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CN ■ I. INTRODUCTION 

cn 

A possible mechanism of Electroweak Symmetry Breaking is provided by strongly-interacting dynamics beyond the 
O , Standard Model (BSM) [H, In this picture, a new strongly-coupled gauge interaction acting at energy scales of 
■ the order of 1 TeV is conjectured. This interaction embeds the Standard Model gauge group SU(2)l U{1)y and 
. . \ contains fermionic degrees of freedom (different from the Standard Model fermions) that do not interact directly with 
K>" • the latter. The chiral symmetry of the BSM interaction breaks spontaneously at the scale of 1 TeV. This breaking 
yTa ] provides mass to the Z and bosons of the Standard Model. Historically, this framework is known as Technicolor; 

. the new bosons are referred to as technibosons and the new fermions as technifermions. In order to give mass to 
d ' the Standard Model fermions, another new gauge interaction acting at higher energy scales (Extended Technicolor) is 
introduced. 

In the original proposals the strongly interacting BSM dynamics was obtained by rescaling QCD. The ratio between 
the scale of the technicolor model and the QCD scale can then be used to determine the mass spectrum of the BSM 
theory. However, this scenario proves to be inadequate to provide a mechanism of mass generation for fermions 
without running into problems with flavor changing neutral currents. A more refined framework that could avoid 
those problems is Walking Technicolor ^3(-t5)]. 

Walking theories are realized as deformations of theories with an infrared (IR) fixed point, i.e. a point in which the 
/3 functions for the couplings of the theory vanish [H, [7| . The role played in this scenario by the fermion representation 
has been emphasized in Refs. In particular, for theories in the two-index symmetric and adjoint representations, 

an IR fixed point can be reached at smaller values of the number of fermion flavors n/ than for theories involving 
fundamental fermions. 
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Like any other BSM framework, Technicolor has to confront the stringent experimental bounds for new physics 
summarized in the S and T parameters [loj . Recent reviews of the phenomenological aspects of Technicolor theories 
can be found in Refs. [Ill4l5l |. It is currently an open question whether a theoretically consistent framework for 
Electroweak Symmetry Breaking can be drawn from those ideas. In particular, one would like to explore from first 
principles whether SU(7V) gauge theories with Uf fermion flavors in the fundamental or in a two-index representation 
can provide a viable walking scenario for some values of N and nj. Such theories would be natural candidates as 
models of strongly interacting BSM dynamics, which eventually will be tested at the LHC. 

Ultimately, the issue of determining the features of a SU(-/V) gauge theory coupled with Uf fermion flavors trans- 
forming according to some representation TZ of the gauge group is of a nonperturbative nature, and as such it can be 
studied in the framework of lattice gauge theories (see e.g. for a complementary approach based on AdS/CFT 
techniques). Following the work of Ref. [17|, other lattice studies have focused on theories with fermions in two-index 
representations conjectured to be relevant as models of strongly interactin g B SM dynamics: preliminary results have 
appeared for SU(2) with two fermions in the adjoint representation |18l422l |. and SU(3) with two fermions in the 
symmetric representation . These studies found a mass spectrum characterized by the degeneracy of the pscu- 
doscalar and the vector meson in a wide range of fermion masses. Complementary investigations of the running of the 
coupling (23 - |29| . and of the exponents that govern the scaling towards the massless limit [30l - l33| found preliminary 
indication for the existence of an infrared fixed point. Note that while in SU(2) with adjoint fermions there seems 
to be consensus that the theory has an infrared conformal fixed point, for SU(3) with sextet fermions the study of 
Ref. [S^ supports a QCD-like scenario, while Ref. [H, favors a walking scenario. A closely related line of research 
is the lattice determination of the critical number of flavors for the onset of the conformal window in the SU(3) gauge 
theory with fermions in the fundamental representation [37l - l48j . Ref. [49| provides a perturbative determination 
of renormalization constants relating lattice and continuum observables in SU(iV) gauge theories with fermions in 
two-index representations. A recent account of the activity in the fleld is provided in Refs. [H, [10, HH. For current 
numerical studies of conformal gauge theories, it proves to be helpful to have analytical estimates of the extent of the 
conformal window; recent works on this subject are reported in Refs. (52| - [55j . 

Numerical simulations of the spectrum of candidate theories of Electroweak Symmetry Breaking beyond the Stan- 
dard Model have focused almost exclusively on the meson spectrum (and in particular, on the states that in QCD 
are the lowest-lying particles of the meson isovector spectrum, namely the pseudoscalar and the vector mesons). Re- 
cently, investigating the case of a SU(2) gauge theory with two fermion flavors in the adjoint representation (which 
is commonly referred to as Minimal Walking Technicolor), we have pointed out in Ref. [56] that a clean signature 
of conformality in the chiral limit can be obtained by comparing mesonic and gluonic observables. In particular, a 
conformal gauge theory broken with a small fermion mass term displays the phenomenon of hyperscaling and locking, 
i.e. all the ratios of spectral quantities are independent of the fermion mass if the latter is sufficiently small. This 
paper has the twofold motivation of discussing more extensively the general expectations for the spectrum when a 
small mass term breaks explicitly conformal invariance, and of presenting the details of our analysis of the data in 
the gluonic sector leading to the conclusions of Ref. [s^ about the likely existence of an infrared flxed point, having 
increased the statistics at some values of the lattice parameters. This work complements the investigation reported 
in Ref. [s^l, where our results for mesonic observables were discussed. We shall use the evidence found for the theory 
to have an infrared fixed point to perform a scaling analysis of our observables as a function of the fermion mass 
and provide an estimate for the anomalous dimension of the condensate, whose value has relevant phenomenological 
implications. 

As in Refs. [s^ [s^], the study reported here is at fixed lattice spacing. One key issue that should be carefully 
discussed is whether our results are relevant for the continuum physics. In general, the program of extracting the 
values of observables in the continuum from lattice simulations of BSM models is still at an early stage; in practice, 
numerical results are obtained for volumes and lattice couplings that are argued to be a good approximations of 
the continuum system. Assessing the reliability of lattice simulations for continuum physics requires then a detailed 
knowledge of the phase structure of the lattice theory. More in detail, it is easy to prove analytically that deep in 
the strong coupling phase a SU(iV) lattice pure gauge theory is always confined. This feature survives when fermions 
with sufficiently high mass are added to the action. The lattice strong coupling, also known as the bulk phase, is 
separated from the continuum phase by either a phase transition or a smooth crossover (for a study of the strong 
coupling regime with fermions in the two-index representation, see [58]). Note that the latter phase may or may not 
be confining. In order to obtain a reliable continuum extrapolation, only points for which the system is in the same 
phase as the continuum theory must be considered. Hence, one of the preliminary tasks of lattice simulations is to 
identify the exact extent of the bulk phase. This program has been carried out in Refs. [18n-20j. which have shown that 
the bulk phase roughly corresponds to lattice couplings /3 < 2.0. However, staying clear from the bulk phase could 
not be a sufficient condition for getting relevant results for the real-world physics: another aspect that needs to be 
considered is the physical size of the volume, which should be such that analytical predictions for finite size corrections 
could be reliably used to extract information from the data. Although at first sight this issue could seem more under 



3 



control than the one related to lattice artefacts, this is a prejudice modeled after our understanding of QCD. If the 
physics of our system is conformal, QCD does not provide a reliable guidance for analyzing the numerical results. In 
fact, recent analytical and numerical arguments have shown that, even in the continuum, the theory in a finite box 
is characterized by a non-trivial phase structure in terms of the lattice volume and of the fermion mass [59l - l6l| . The 
finite volume phases are characterized by different behaviors of the order parameter related to confinement, the trace 
of the Polyakov loop. In order to check for possible lattice artefacts, we measured the Polyakov loop wrapping around 
the spatial directions. We find that the bare fermion mass can significantly affect the phase structure at fixed lattice 
size. In particular, at fixed volume, when lowering the bare fermion mass the system goes from the infinite-volume 
confined phase (as shown by the presence of a non-zero string tension) to a spatially deconfined regime. We find that 
gluonic observables are strongly affected by this change of regime, while mesonic quantities behave smoothly across it. 
This could be an effect of the separation at finite mass of the confinement and the chiral symmetry breaking scales in 
theories with adjoint fermions (observed in Refs. [62l. l63j). which would be an expected feature for a candidate model 
of Technicolor that has a large distance dynamics different from the QCD one. 

This work is organized as follows. Sect. |ll] discusses the hyperscaling and locking phenomena, which arises when 
an infrared conformal gauge theory is deformed with a small mass term. In Sect. IIIll we define the discretized theory 
and set the notations. Results for the phase structure, the string tension extracted via Polyakov loop correlators, the 
string tension extracted from expectation values of Wilson loops and glueball masses are reported in Sect. IIVIIVIII 
A comparison between the dynamical and the quenched simulations performed keeping the string tension and the 
pseudoscalar mass fixed at the values dictated by the dynamical theory is then provided, and from this comparison 
hyperscaling and locking are shown to take place in the model studied in this work (Sect. IVIIII) . Sect. IIXI illustrates 
our scaling analysis aimed to determine the chiral condensate anomalous dimension 7^,. Finally, Sec. 1x1 reports our 
conclusions and possible future directions of our work. 

II. MASS-DEFORMED INFRARED-CONFORMAL GAUGE THEORIES 

For technical reasons which depend on the specific fermion discretization, lattice simulations can only be performed 
with a non- vanishing mass term for the fermions. In particular Wilson fermions break chiral symmetry explicitly even 
for vanishing bare mass, so that the massless limit is only obtained by fine-tuning the parameters in the Lagrangian. 
Moreover, it is impossible to simulate at arbitrarily small masses if the lattice spacing and the volume are kept 
constant, since small eigenvalues of the Dirac operator are generated, the simulation algorithm becomes unstable and 
unphysical phases can appear. The extrapolation from a region of small enough masses (but still in a safe region of 
parameters) to the chiral limit can be performed only under the guidance of an analytical picture. 

For QCD-like theories, chiral perturbation theory (in the infinite volume, in the epsilon and delta regimes) allows 
to extrapolate physical quantities from a region of small enough masses to the chiral limit. 

The natural question is: what should we expect if we deform an IR-conformal theory with a small mass term, and 
how do we recover the chiral limit? Hence, before illustrating the details of our simulations, we set the frame for the 
picture in the latter case. 

A. Renormalization Group analysis 

A gauge theory with massless fermions (in the continuum) depends on a single parameter, the running coupling 
g{lJLo) at some reference scale /io, or alternatively the RG-invariant parameter A. This is valid for both confining and 
IR-conformal theories. In confining theories the particle masses (except the Goldstone bosons) are proportional to the 
parameter A in the chiral limit. In the case of IR-conformal theories, where the spectrum is made of unparticles, A 
is not associated to particle masses, but sets the energy scale at which the cross-over between the asymptotically free 
and conformal regimes occurs. An explicit definition of A is not relevant for our discussion and will then be omitted. 

When the IR-conformal theory is deformed by a mass term for the fermions, a particle spectrum with a mass gap is 
expected to be generated. The theory depends now on one more parameter, the running mass m{p,Q) at the reference 
scale ^0; alternatively, an RG-invariant parameter M can be suitably defined. Close enough to the chiral limit (in 
the scaling region), the particle masses are expected to be independent of A. We will see that under a regularity 
hypothesis, those masses are expected to be proportional to AI. This result is standard in the statistical-mechanics 
analysis of second-order phase transitions, but it will be presented here using the language of quantum field theory. 

The running of the renormalized mass is computed by solving the RG equation (in a mass-independent renormal- 
ization scheme): 



= -7(.9(m))"^(m) 



(1) 
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which yields: 



i(^) = m(^o) cxp -( - / — — -dz ^ = Z,„(//,/io, A)m(/xo) . (2) 



9(A'o) 



As we are going to show, the function Z„i(/i, /io, A) can be rewritten in a more convenient form. The theory we are 
interested in is asymptoticaUy free in the UV. The (3 and 7 functions close to the UV fixed point are: 

g^O : Pig) ~ -/Jqs' , (3) 

7(5) - 705' , (4) 

where the lowest order coefficients come from a one-loop computation (Tr and C2{R) are the generator normalization 
and the Casimir of the fermionic representation): 

/3o=(^(^A^jT,n.) , (5) 

^From now on, we will be interested only in the IR-confornial scenario. Close to the IR fixed point we assume a 
regular behavior for the RG functions: 

5-^.9* ■ Pig) - P*{9 - 9*) , (7) 

7(5) - 7* , (8) 

where and 7*, which are scheme-independent quantities, are in general not accessible by a perturbative expansion. 

Integrating the RG equation iidg/d/j, — Pig) close to the fixed points, the asymptotic running-coupling behavior is 
derived: 

H^oo : gifi) ~ ■ ^, ... , (9) 
2polog(/i/A) 

^^0 : gi^i) ^ g, - Ag i^^Y' . (10) 

We separate now the singular behaviors close to the fixed points in the multiplicative renormalization function of 
the mass: 



- / 

"'9(^0) 



7(f) _ 7* _^ Ja_ 
Piz) P^iz-g^) Poz 



dz } X 



^^^^ 7* ^ I j r^''^ 70 



X exp <^ - / in ^ exp <^ / —dz 

Zynin/K) 



Zm(^o/A) 

where the function 



(11) 
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ZMK) = [5, - exp J /'* f ^ - + ^) ^^4 (12) 

[^(m) VP(^) PAz-9*) PqzJ J 

is defined in such a way that the integral in the exponential is finite both for /i — > and /i — > 00. 

An RG-invariant fermionic mass M can be defined by means of the condition m(M) = M . Plugging Eq. (ITT]) in 
Eq. ©, and choosing = M we get the relationship: 

ZMA)-'mifi) = Z„,(M/A)-iM . (13) 

If the RG-invariant mass M is much larger than A, the following asymptotic behavior can be easily shown to hold 
by using the previous relationship: 



mip) = AooZ„,i^i/A)M 



1 ^ 



(14) 
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At fixed energy scale /z, the running mass diverges as M goes to infinity. The fermions decouple and the theory is 
effectively described by a pure Yang-Mills with a scale Aym — A. At leading order in A/M, the parameter M coincides 
with the quark pole mass. In fact, if S{p) is the quark propagator in a fixed gauge, the perturbative expansion yields: 



sip)- 



71=1 



(15) 



the pole mass m is defined in such a way that the quark propagator has a pole for — — . The pole mass is RG 
invariant, therefore it can be computed for an arbitrary value of /i. It is convenient to choose /i = M: 



2 

\ ^ „ / rn 



^^ + l.^"(l72'l)5'"W- (16) 

n=l ^ ^ 

At large masses M 3> A, the terms in the sum are suppressed since the running coupling goes to zero, and m ~ M . 
In this regime, the meson masses are just twice the quark pole mass, while the glueball masses are the same as in the 
pure Yang-Mills theory: 

M^es = 2M ; (17) 
-^:^giuc — ^giucA . (18) 
On the other hand, the chiral limit is reached for values of M much smaller than A. In this case, Eq. (IT51) becomes: 

m(^) = AoZ(^/A)A-^-Afi+T* , (19) 

producing the power law that is characteristic of the IR fixed point deformed with a small fermionic mass. 

Consider now a physical mass Mx in a channel X (it can be the mass of a particle or other physical quantities like 
the square root of the string tension). As every observable, this will be a function of the renormalized coupling g{^)^ 
the mass m{fi), and the subtraction scale /i. However a physical quantity must be RG invariant: 

Mx[t^,g{ti),m{ti)]= Mx (20) 

for every value of /i. The RG equation for Mx has a simple solution in terms of the RG-invariant quantities A and 
M: 

Mx=M Fx{M/k) , (21) 

where Fx is a generic function of the ratio M/A. In particular, if Fx{x) = a/x, we get Mx = a A which is an 
RG-invariant quantity, but does not vanish in the chiral limit. 

The hyperscaling hypothesis, which is assumed in the standard discussion of second-order phase transitions (see e.g. 
Ref. [S^l), asserts the regularity of masses (or correlation lengths in the language of statistical mechanics) with respect 
to the irrelevant couplings. Consider Eq. ([^0]) for /i = M <C A: 

Mx ^Mx [M, - Ag{M/Af' , M] = 
=Mx[l,9.- Ag{M/Af%l] A/~ 

^Mx[l,9.A]M = AxM, (22) 

where we used dimensional analysis for the second line, and regularity with respect to g in the last one. Under the 
hyperscaling hypothesis, RG-invariant IR quantities depend only on M (and not on A) close enough to the chiral 
limit. The hyperscaling hypothesis constraints the Fx function defined in Eq. (I^Tt to be regular in the chiral limit: 

\iuiFx{x)^Ax. (23) 

Since a mass gap is expected to be generated at nonzero values of M, Ax must be different from zero. 
Combining Eqs. ((22|) and (|T9l) we get the power law for physical masses close to the chiral limit: 

Mx = AxM = Ax[AoZ(/^/A)]"TT^ATT^m(^)TT^ . (24) 

We remind that this expression is valid for every value of /i as long as Af <C A. In particular, the independence of 
Mx of A is manifest at values <C A: 

Mx = Ax^^'^^m{^l)^^ . (25) 

If we interpret the RG in the Wilsonian sense and choose fi — a^^ to be the cutoff, Eq. ([M)) yields the power law 
dependence of physical masses on the bare quark mass aAlx o: (amo) . 
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FIG. 1: Sketches of the spectrum of a mass- deformed IR-conformal theory (square root of the string tension, 0^"*" and 2^"*" 
glueballs, pseudoscalar and vector isovector mesons). In the left plot, the locking sets up at an intermediate value of the fermion 
mass, where dynamical fermion effects account for the physics of the system, but the pseudoscalar is not much lighter than the 
other particles in the spectrum. In the right plot, the locking sets up at a high value of the fermion mass, where the heavy quark 
effective theory provides a good description of the relevant degrees of freedoms. This case is realized close to the Banks-Zacks 
point, but is possible in principle also if a strongly coupled IR fixed point is present. 



B. Scaling region and locking scale 



Under the hyperscaling hypothesis, the function Fx defined in Eq. (I21[) is expected to approach a nonzero value 
Ax in the chiral limit. We can define the scahng region for a given channel X as the range of a; = Af/A around a; = 0, 
where the function Fx{x) deviates from its asymptotic behavior by a small relative amount e: 



Fx{x)-Ax 



Ay 



< e . 



(26) 



In the scaling region, the mass Mx obeys the power law ([M)) as a function of the running mass up to corrections 
of order e. The extension of the scaling region will depend on the size of the discarded subleading contributions to 
formula ([M]) in the chosen channel. 

Consider now the square root of the fundamental string tension Mg- = ^fa (which is well defined for dynamical 
fermions in the adjoint representation) and the lightest isovector meson (which is always the pseudoscalar one), with 
mass Afpg. A finite value x — x exists, below which both these channels are in the scaling region. This means that 
below the mass Miock = a;A, the corrections to the hyperscaling behavior of M„ and Mps masses are relatively smaller 
than e. Also the ratio M-ps/^a for every fermionic mass below Miock will be very similar to its asymptotic value 
Avs/Aa-. 



Mps Aps 



Ma Aa 



<0{e) 



(27) 



The dynamics is dramatically different below and above the mass Miock- In the large-mass region, M ^ A, the 
gluonic and mesonic masses are parametrically independent. All the gluonic masses are proportional to A, while all 
the mesonic masses are equal to 2M: 



Mps = 2M , 
Ma = BaA . 



(28) 
(29) 



The ratio Mps/Ma- goes to infinity in the large- mass limit. For masses below Afiock the two masses Afps and M^ 
enter the scaling region, become both independent of A and proportional to M. The ratio Afps/A/o- is locked to its 
asymptotic value Aps/Aa-- We will refer to Afiock as the locking mass. 
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The behavior of the masses in between the large-mass and scahng regions and the actual value of x depend on the 
details of the dynamics. However if the dynamics is such that the locking occurs at a value x = Miock/A ^ 1, then 
both Ala- and Afps at the locking scale are still approximately the same as in the large-mass region: 

M,{M = Afiock) ^^A/iock _ , _ 

i>(T — — — ^crX , [OU) 

Aps ^ 2 , (31) 
and the ratio Mps/M^ is locked at a very large value: 

Aps 2x 



> 1 . (32) 



Mesons are much heavier than the square root of the string tension for every value of M . Choosing an intermediate 
energy scale E such that <^ E <^ Mps, the effective theory describing the gluonic degrees of freedom at energies 
below E is a pure Yang-Mills plus power-suppressed corrections coming from the propagation of heavy quarks in the 
loops. In order to write the effective Lagrangian in this regime, we need all the gauge-invariant scalar operators of 
dimension 6 that are invariant under parity, and charge conjugation. These can be written as linear combinations of 
the following independent operators (a similar analysis on the lattice was carried on in Ref. (65|): 

Si = ^ tr ( J^.p J'^"") , (33) 

^2 = 5] tr ( J'^^p J/^) , (34) 

S3 = ^ tr ( J^,prn , (35) 
where Jpi,p = d^^F^p — i[Afj_, F^p\. Thus the effective Lagrangian can be written as: 

'Ceff = -^tr(Fp.Fn+ E + (36) 

^ 1=1,2,3 

The scale Aym of this low-energy pure Yang-Mills is in general a function of A and M and can be computed by 
matching the square root of the string tension of the low-energy effective theory with the same quantity computed in 
the dynamical theory: 



YM 



1 + 



Aym^ ^ 
M 



^M, = MF,{M/k) , (37) 



which implies that trivially Aym — A for M ^ A, while for M < M\ock then 

Aym ^ ^ — ^ — ■ (38) 

In the scaling region the scale Aym of the low-energy pure Yang-Mills slides with the RG-invariant fermionic mass 
M. 

A comment is mandatory at this point. At fixed value of the fermionic mass, the low-lying spectrum of a mass- 
deformed IR-conformal theory with a: ^ 1 can not be distinguished by the low-lying spectrum of a confining theory 
with heavy quarks, since they both are described by the same effective Lagrangian p6p . However in a genuine heavy- 
quark phase the low-energy spectrum is almost independent of the mass M, while the sliding of the low-energy scale 
described in Eq. psp and (equivalently) the locking of the gluonic spectrum to the mass M is ultimately a very clean 
signature of IR-conformality. 

Summarizing: 

• We define the locking mass Miock as the mass below which both the lowest isovector meson and the string 
tension are approximately in the chiral scaling region. 

• The value of x = Mjock/A is determined by the detailed dynamics of the theory. If a; ^ 1 then the mesons 
are always much heavier than the square root of the string tension. The low-energy effective theory is a pure 
Yang-Mills plus small corrections, with a scale Aym which depends on both A and M. For M > Miock ^ A 
then the fermions completely decouple and Aym — A, while for M < Miock the only effect of the fermions in 
the dynamical theory is to make the low-energy scale slide with the fermionic mass Aym — M/x. 
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• The case where x 3> 1 is reahzed if the fixed point is perturbative [66l |67|. In fact, in this case 



(39) 



The described scenarios are ihustrated in the sketches in Fig. [TJ 



III. THE LATTICE MODEL 



Consider a four-dimensional Euchdean torus L^x L^, where Lt and Ls are the lengths respectively of the temporal 
and spatial directions. The space-time is discretized by introducing a lattice with spacing a, and with Nt = Lt/a and 
Ng — Ls/a sites respectively in the temporal and spatial directions. Lattice sites are identified by four-coordinate 
dimensionful vectors x = (xq, xi, X2, xa). Therefore xo/a is an integer number running from zero to Nt — 1, and 
Xi/a with i — 1,2,3 are integer numbers running from zero to Ns — 1. In some cases it is useful to separate the 
temporal coordinate from the spatial vector; we write x = (t, r), and r is the modulus of r. Lattice directions are 
indicated with a Greek symbol and run from to 3. The temporal direction is chosen as the zero-th direction. We 
will use the same Greek symbol both for the direction index and for the vector of length a along the axis direction 
(the meaning of the symbol will be always clear from the context). The spatial directions are closed with periodic 
boundary conditions (PBC) for all fields, while the boundary conditions in the temporal direction are periodic for 
gauge fields and antiperiodic (ABC) for fermion fields. 

The action of a SU(Af) gauge theory with fermions can be decomposed as 



where Sg is the discretized Yang-Mills action and 5*/ is the fcrmionic contribution. Various choices for the lattice 
action are possible, differing from each other by corrections that vanish in the continuum limit. At finite lattice 
spacing different choices are differently affected by lattice artefacts. In particular, as the lattice spacing is increased, 
a transition to a phase not connected with the continuum (the bulk phase) takes place. A careful exploration of the 
phases of the system on a lattice as a function of the lattice parameters is then mandatory. 
For the gauge part, we use the Wilson action: 



where V^v{x) is the parallel transport of the link variable U (a;, /i) e SU(A^) along the elementary square of the lattice 
identified by the point x and the pair of directions (^, v). (3 is related to the bare coupling g'^hy (3 = 2N/gQ. The value 
of the coupling determines the physical value of the ultraviolet cut-off, the lattice spacing a. Note that independently 
of the fermion representation, the link variables are in the fundamental representation of SU(A^). 

The fermion part of the action for a spinorial field '4'{x) defined on sites x and transforming in the representation 
R can be written as 



(40) 




(41) 




(42) 



k=l 



where Dm is the Dirac operator, in the Wilson discretization: 



(43) 



where U are the link variables in the representation R, and mo is the bare mass. 
The functional integral is given by 




(44) 



and the vacuum expectation value of an operator 0{U, tpi V') by 




(45) 
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where once again it is possible to integrate over the fermion fields and obtain an expression that involves only an 
integral over the link variables. For further details on the lattice formulation, we refer to Refs. [isllSTj. 

We performed numerical simulations for SU{2) gauge theory with Uf — 2 Wilson fermions in the adjoint repre- 
sentation at fixed value oi (3 = 2.25, different values of the bare mass, and different lattices, the smallest one being 
a 16 X S'^ lattice, and the largest one being 32 x IQ"^. We used the RHMC algorithm [6^ as implemented in the 
HiRep code, which is described and benchmarked in detail in Ref. fis'l. The full list of the parameters we used in 
our simulations can be found in Tables Hl lIIIl In this work we are mainly interested in gluonic observables. For every 
choice of the parameters we compute: 

• the traced Polyakov loops in every direction, in order to identify the regime of the theory, as described in 
Sect.HV] 

• the string tension by means of correlators of spatial and temporal Polyakov loops, as described in Sect. |Vl 

• the static force and potential, as described in Sect. IVIl 

• the glueball masses, as described in Sect. I VIII 

We follow the convention that lattice observables are dimensionful, with the same dimension of the corresponding 
continuum observable. Of course we can measure only dimensionless ratios. For instance given a mass m, only the 
dimensionless quantity am can be extracted from lattice simulations. Determining a(/3), and ultimately the physical 
value m, requires to set the physical scale using an appropriate observable. We will not perform the step of reinstating 
physical units, but shall leave a (which is fixed in our case) as a parameter. However, the reader must bear in mind 
that knowing the value of a in terms of the quantities entering the dynamics in the continuum is important in order 
to confidently assess the relevance of a lattice simulation for continuum physics. For this investigation, following the 
detailed exploration of the phase structure of the theory performed in Ref. [l^j, we argue that large discretization 
artefacs are ruled out and we postpone to future studies a systematic investigation of these effects. 



lattice 


V 


—amo 


spatial center 
Nconf realization 


string tension from Polyakov correlators 

1/2 1/2 

aa^ aus as = Ut 


SO 


16 X 8^ 


-0.5 


8000 


S 






SI 


16 X 8^ 


-0.25 


8000 


s 


0.4085(57) 0.393(11) 


yes 


S2 


16 X 8^ 


-0 


8000 


s 


0.3998(57) 0.388(11) 


yes 


S3 


16 X 8^ 


0.25 


8000 


s 


0.328(23) 0.358(12) 


yes 


S4 


16 X 8^ 


0.5 


8000 


s 


0.3576(46) 0.347(11) 


yes 


S5 


16 X 8^ 


0.75 


8000 


s 


0.282(13) 0.2784(75) 


yes 


S6 


16 X 8^ 


0.90 


8000 


s 


0.227(11) 0.2452(67) 


yes 


AO 


16 X 8^ 


0.95 


8000 


s 


0.1974(83) 0.2218(35) 


no 


Al 


16 X 8^ 


0.975 


8000 


? 


0.2066(97) 0.2094(49) 


yes 


A2 


16 X 8^ 


1 


8000 


? 


0.1960(85) 0.2252(62) 


no 


A3 


16 X 8^ 


1.025 


8000 


A 


0.1689(44) 0.2109(46) 


no 


A4 


16 X 8^ 


1.05 


8000 


A 


0.1679(47) 0.2074(38) 


no 


A5 


16 X 8^ 


1.075 


6400 


A 


0.1629(27) 0.20680(94) 


no 


A6 


16 X 8^ 


1.1 


6400 


A 


0.1553(28) 0.20443(82) 


no 


A7 


16 X 8^ 


1.125 






0.1462(26) 0.20423(87) 


no 


A8 


16 X 8^ 


1.15 






0.1368(20) 0.20402(67) 


no 


A9 


16 X 8^ 


1.175 






0.1310(19) 0.2067(11) 


no 



TABLE I: Results for the Polyakov loop distribution and the string tension, for all the simulations on the 16 x 8^ lattice. 
We analyzed the Polyakov loop distribution by using the number of configurations quoted in the 4th column (for the three 
lowest masses we could not safely discard the thermalization) . In the 5th column, 'S' indicates that the distribution has a 
single maximum in zero, 'A' indicates that a double peak or an asymmetric peak is visible, '?' indicates that we cannot 
clearly distinguish between the two cases. The temporal (6th column) and spatial (7th column) string tensions computed from 
correlators of Polyakov loops are quoted. In the 8th column, 'yes' indicates that the data for the temporal and spatial string 
tensions have an overlap at la. 
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lattice 


V 


—amo 


spatial center 
Nconf realization 


string tension from Polyakov correlators 

1/2 1/2 

acTj aas as = o"t 


RD 

JDU 


94- V 19'' 


n 


9501 


S 


0.225(11) 0.2181(23) 


yes 


Bl 


24 X 12'' 


1 


/9oi 


b 


n i'70c;^c;'7\ A iooo//in\ 
U.L I00[0 1 ) U.i802(4U) 


yes 


B2 


24 X 12'' 


1 0'^ 


5620 


"? 


0.1491(73) 0.1597(24) 


yes 


B3 


24 X 12'' 


1.075 


4987 


A 


0.1398(63) 0.1461(20) 


yes 


RzL 


9/1 V 19'' 


1 1 
1 . 1 


4194 


A 


0.1205(56) 0.1427(21) 


no 


B5 


24 X 12^ 


1.125 


4001 


A 


0.1130(54) 0.1338(5) 


no 


B6 


24 X 12^ 


1.15 


1500 


A 


0.0829(26) 0.1337(9) 


no 


B7 


24 X 12^ 


1.175 


5001 


A 


0.0918(36) 0.1333(6) 


no 


B8 


24 X 12^ 


1.18 


4490 


A 


0.0944(36) 0.1333(11) 


no 


B9 


24 X 12^ 


1.185 


4335 


A 


0.0834(29) 0.1425(47) 


no 


BIO 


24 X 12^ 


1.19 


4336 


A 


0.0851(34) 0.1380(37) 


no 




TABLE II: As Table H] for the simulations on the 24 x 12^ 


lattice. 


lattice 


V 


—amo 


spatial center 
Nconf realization 


string tension from Polyakov correlators 

1/2 1/2 

acTj aas as ~ at 


CO 


32 X 16^ 


1.15 


6145 


? 


0.0790(34) 0.1029(19) 


no 


CI 


32 X 16^ 


1.175 


1871 


A 


0.0966(78) 0.10057(65) 


yes 


C2 


32 X 16^ 


1.18 


1500 


A 


0.0648(33) 0.1086(13) 


no 


C3 


32 X 16^ 


1.185 


1419 


A 


0.0612(36) 0.09953(13) 


no 


C4 


32 X 16^ 


1.19 


1609 


A 


0.0703(40) 0.1021(13) 


no 



TABLE III: As Table U for the simulations on the 32 x 16^ lattice. 



IV. CENTER SYMMETRIES AND POLYAKOV LOOPS 

For a SU(A^) gauge theory with fermions in the adjoint representation, the action has a invariance, where each of 
the Zjv factors is associated with one direction of the system. An observable that is not invariant under the Zjv factor 
of the symmetry group associated to the direction p is the Polyakov loop operator in the fundamental representation 
wrapping around p: 

Pp^T. ti Uix + np,p) , (46) 

where the primed sum runs over all points of the lattice slice at Xp = 0. Consider the system in a Euclidean manifold 
E'^ X S"^, in which the direction p is compactified on the circle and the other three directions extend to infinity. 
At given radius of the S^, the Ztv symmetry associated to p might either be a symmetry of the system or could be 
spontaneously broken. For simplicity, let us take p to be the temporal direction; the inverse radius of the has 
then the natural interpretation of the temperature of the system. In this case, if the system is symmetric, the trace 
of Pp is equal to zero and the system is confined at the given temperature. Conversely, a non-zero (trPp) signals the 
breaking of the Zjv symmetry associated with the direction p and the system is deconfined. Analogously, if p is a 
spatial dimension, we call the broken phase spatial deconfinement. 

Even if the theory is conformal in the chiral limit, a small mass drives it off the attraction basin of the IR fixed 
point, generating a mass gap. As for the pure gauge system, at finite temperature the theory is expected to undergo a 
deconfinement transition as the temporal direction is reduced down to a critical value. Due to the periodic boundary 
conditions for fermions in space, the reduction of a single spatial direction is more subtle. Recent numerical simula- 
tions 60] have shown that in a SU(3) theory with two staggered Dirac adjoint fermions the center symmetry (which 
is intact at infinite volume) is first broken and then restored again as the radius of a compactified spatial direction 
is shrunk. The symmetry restoring transition happening when the radius is shrunk from values that set the system 
in the broken phase is expected from the one- loop perturbative calculation of Ref. [69]. The critical radii for the 
symmetry breaking and symmetry restoring transitions depend in general on the mass. Although the different dis- 
cretization choice and the different gauge group might affect the phase structure, it is possible that a similar behavior 
characterizes SU(2) with two adjoint Dirac flavors of Wilson fermions. 
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When the system is in a compact domain, non-trivial phases can be observed only at large N. In this case, a rich 
phase structure can exist, since more than one Z^r can be broken at the same time. For the pure gauge theory at 
large N, as the hypercubic volume is reduced from large values, a cascade of phase transitions takes place ^7^. Each 
phase can be characterized by the number of the Zjv factors that are spontaneously broken to the trivial element of 
the group. The case of a large- SU(Ar) gauge theory with adjoint fermions on a compact domain has an even richer 
phase structure [7l|. 

If the three spatial directions have a finite extension and N is finite, no phase transition can occur at any size of the 
system. However one can still investigate whether the distribution of the Polyakov loop in a given spatial direction 
has a single peak in zero (S, symmetric phase) or two separate peaks symmetric around zero (A, asymmetric phase), 
keeping in mind that at infinite volume the Polyakov loop distribution should display a single peak. The S-phase is 
the finite volume and finite N equivalent of the thermodynamic phase in which the system is symmetric under the 
Ztv symmetry related to the direction wrapped by the Polyakov loop, while the A-phase is the finite volume and finite 
N equivalent of the broken phase. We stress once again that at finite N and on a finite volume there are no distinct 
phases, but only different regimes. The terminology here is used only for convenience. 

We measure the distributions of the Polyakov loops in all the directions. In all the lattices we consider, we always find 
a temporal S-phase, which means that we are correctly simulating the confined thermal phase of the infinite volume 
system. In lattices with geometry Nt x Nf, we are not interested in separating the equivalent spatial directions. Hence 
we will say that the system is in a spatial S-phase if all the spatial directions show a single peak distribution centered 
in zero. 

For the 16 x 8^ lattice, we find the spatial A-phase below amo « —0.975. The crossover from the S-phase to the 
A-phase is very smooth. A summary of the realization of the spatial center in the 16 x 8^ can be found in Table ID 
The two maxima of the spatial distributions are better defined on the 24 x 12^ lattice, for which the spatial A-phase 
is found below auiQ = —1.05 (see Table HIl for a summary). For the 32 x 16^ lattice, all the simulated masses show 
distinct maxima in the distribution of the spatial Polyakov loops, except at the higher mass where we find a broad 
distribution with a flat region in the middle (see Table iHll for a summary). We show some distributions for the 
24 X 12^ and 32 x 16'^ lattices in Fig.O We also simulated a 24 x 24 x 12^ lattice at amo = —1.125. In this case the 
large spatial direction is in the S-phase, while the smaller ones are in the A-phase. In Fig. |3]we plot the histories and 
the distributions for all the Polyakov loops for the 24 x 24 x 12^ lattice. 

In Fig. m we also report the absolute value of the Polyakov loop on the 16 x 8^ and 24 x 12'^ lattices. As the 
mass is decreased, this quantity undergoes a crossover from a lower value to a higher one. This crossover moves to a 
noticeably lower mass when the lattice size is increasing, indicating that the A-phase disappears in the large- volume 
limit. 

The crossover from the S-phase to the A-phase as the spatial volume is decreased is a clear indication that at least 
our S-phase is not in the femtoworld regime, which is in the nearby of the zero-volume limit, and which is a possible 
source of large systematic errors in a lattice simulation [l^ . Our data are actually consistent with the picture that 
the S-phase is connected with the infinite volume limit. However, as we shall show in detail, we still have large 
finite- volume effects for several of the measured observables. 



V. STRING TENSION FROM CORRELATORS OF POLYAKOV LOOPS 

For a SU(iV) gauge theory in the confined phase, a static quark- antiquark pair in the fundamental representation 
at large separation R is bound by the potential 

V{R) = aR , (47) 

where the string tension a is the dynamically generated scale of the system, a is the string tension in an effective string 
theory describing the low energy dynamics of confining flux tubes connecting the quark and the antiquark. Contrary 
to the adjoint string tension (and to the fundamental string tension in QCD), the chromoelectric field between two 
fundamental sources in a gauge theory with adjoint matter is not screened. Hence, the asymptotic fundamental string 
tension is a well-defined quantity. 

It is easy to prove analytically on the lattice that any gauge theory has a non-zero string tension at strong coupling 
and large fermion masses. The relevant question for the system under study is whether a region in bare parameter 
space exists, which is analytically connected with the continuum limit, and where the string tension is zero in the 
massless limit, as it should be if the theory is conformal. 

The string tension can be extracted from correlators of Polyakov loops. In particular, consider the plane defined by 
xo = and a fixed transverse coordinate in one arbitrary spatial direction (e.g. we can consider the case of constant 
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24x 12^ 
am. = -0.950 




-0.08 -0.04 0.04 0.08 -0.08 -0.04 0.04 0.08 -0.08 -0.04 0.04 0.08 -0.08 -0.04 0.04 0.08 
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FIG. 2: Some distributions of Polyakov loops in the temporal and spatial directions. The distributions have been symmetrized 
by hand. At fixed volume, the distribution of the spatial Polyakov loop shows a single peak at zero at the higher mass, and 
it develops two peaks at the lower mass in some of the spatial directions. At fixed mass, it shows a single peak on the larger 
lattice, and it develops double peaks on the smaller lattice in some of the spatial directions. 



coordinate Xi in the first direction), and let us define: 

M^i) = ^ ^Tr I II U{ano,r;0) 1 , (48) 



1 1 \ 

]^Tr n U{ano,r;0)\ 

* X2,X3 \no=0 / 



where r is the spatial vector with coordinates {xi,X2,X3). For the vacuum-subtracted correlator of this quantity one 
finds 

oo 

{Po{xiyPo{xi+X^)) - |(Po(a;i))|' = ^ |c„(L,)pe-^"(^')^i , E„ > Vn , (49) 

ri=0 

where the sum runs over all states \n) of the Hamiltonian with non-vanishing overlap c„ — {0\PQ{xi)\n) . At large 
Xi the sum is dominated by the term \co{Lt)\'^e~-^°'^^*'''^^ associated with the exponential with the lowest decay 
rate, Eq^Li), which is the energy of the groundstate in the fundamental string sector. For a confining theory, up to 
subleading corrections one finds [72|| 

Eo{Lt) = aLt - c7r{D - 2)/{6Lt) , (50) 

where c depends on the number of massless fermionic and bosonic modes propagating along the string and D is 
the dimension of the system (_D = 4 in our case). Analogously, if the Polyakov loop wraps a spatial direction, its 
vacuum-subtracted zero momentum correlators in the temporal direction define the so-called spatial string tension. 
For a confined theory at zero temperature, the spatial and the temporal string tensions coincide if the system is the 
S-phase. 

In SU(A^) Yang-Mills theories, the effective theory describing the large-distance dynamics of the confining fiux tube 
is a bosonic string theory [73i - ,7& l . General arguments lead to the following expansion of Eq 77] 
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1000 2000 3000 4000 5000 

FIG. 3: Histories and distributions of Polyakov loops in the temporal (with fermionic antiperiodic boundary conditions, ABC) 
and spatial directions (with fermionic periodic boundary conditions, PBC), for bare mass —1.125 on the 24 x 24 x 12^ lattice. 
The distributions have been symmetrized by hand. 
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FIG. 4: Absolute value of the Polyakov loop on the 16 x 8^ and 24 x 12^ lattices. 

where cg is a term O {{1 / {a^ Lf)) . Note that Eq. ((5T|) is the truncation to second order in l/{aL^) of the ground state 
energy of a bosonic string of the Nambu-Goto type: 



Eo{Lt) = aLtJl- 



Ti{D-2) 



(52) 



This might suggests that the effective string theory describing the large-distance dynamics of the confining flux tube 
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is actually Nambu-Goto. These considerations do not generalize immediately to the case of dynamical fermions in 
the adjoint representation. In fact differently from SU(iV) pure Yang-Mills, in this case it is possible to construct 
explicitly fermionic open-string states. Hence, it could be that the effective theory is not bosonic, in which case 
the firmest result available is given by Eq. (|50p . with c (unknown a priori) counting the zero- modes of the effective 
string. In our analysis, we will assume that the effective theory is bosonic and the string tension will be obtained from 
correlator of Polyakov loops assuming Eq. (j52p . This assumption is justified a posteriori by two crucial observations: 
first, the string tension obtained in this way is in complete agreement with the static potentials and forces computed 
by Wilson loops (Sect. rvT|. for which no effective string theory is assumed; in addition, the low-energy dynamics will 
be found to be an effective Yang-Mills with small corrections (accordingly to the discussion in Sect. IIIBI) for all the 
simulated masses, which shows that the choice of Eq. (|5^ is self-consistent. 

For a theory that confines at zero temperature and undergoes a deconfinement phase transition at some critical 
temperature, the correlator in Eq. (I49[) still decays exponentially with the distance in the deconfined phase, provided 
the sources are screened; however, the inverse of the corresponding energy is now associated with a screening length, 
the Debye screening length. Hence, the exponential decay of Polyakov loop correlators by itself does not imply the 
existence of a string tension. In order to see that the theory is confining, the validity of Eq. ((52|) as Lt is varied needs 
to be proved. Alternatively, one has to show that the static potential (Sect. IVT)) is asymptotically linear and the slope 
of the linear part is related to the ground state mass extracted from Polyakov loop correlators via Eq. (|5^ . 

In general, extracting numerically the string tension from the correlator (j49p proves to be technically hard, since 
the signal-over-noise ratio decays exponentially with the separation. A good degree of success is achieved if the 
zero-momentum Polyakov line is replaced by a fuzzy operator, and a reliable signal can be obtained if a variational 
procedure that involves different fuzzy operators is set up. There are several ways of achieving this; here we follow 
Refs. [tsI. |79|. In practice, a recursive procedure is implemented, which allows to obtain smeared links at step / + 1 
from the fuzzy links at I via the relationship 

(x,Ai)=Proj (x, At) +aS'(') (a:, , (53) 

where S''^'^(a;, /i) is the sum of the four length three non-backtracking lattice paths from x to x + afi (staple) and 
D^^''{x,iJ,) is the sum of the 16 length-five non-backtracking paths with the same start and end points (only the 
directions that are orthogonal to the direction in which correlations are taken enter the sums). The constants a and 
6 are are fixed empirically in such a way that the signal is optimal. Since the (weighted) sum of the paths is not an 
element of the group, to obtain an object that can be interpreted as a fuzzy link this sum needs to be reprojected 
onto SU(A^); this is the meaning of the operator Proj() in Eq. ((55)) . After k steps of smearing, consecutive pairs of 
smeared links going in the same direction can be multiplied to produce blocked links. The combination of smearing 
and blocking yields the link set {U^''\x, ^i)} at blocking level 6, which can be used to compute the fuzzy Polyakov 
loop operator 

1 1 f^'^'' \ 

Po'H^) = J^T.N^' n U^''H^no,r;0) , (54) 

* j/,z ^"0=0 J 

and analogously for the other directions. For sake of definiteness, we discuss the case of Polyakov loops winding in 
time, but similar conclusions hold for Polyakov loops wrapping around the other directions. The element be of the 
correlation matrix C{Xi) is then defined as 

CUXi) = {P|,''\x^)^P|^^\x^+X,)) - {p|^''\x,))^p|,^\x^)) . (55) 

As functions of A"i, diagonal correlators involving the eigenvectors associated with the largest eigenvalues of 
(7~^(0)(7(1) decay as single exponentials with energies Eo,Ei,... already at distances of a few lattice spacings. 
The negligible contamination from excited states eliminates the need to go to large distances to identify the stringy 
state with the lowest energy; at the same time, this procedure provides an estimate of energies of excited states, asso- 
ciated to the single-exponential behavior of diagonal correlators of eigenvectors corresponding to smaller eigenvalues, 
although the reliability of the excited spectrum depends on how large the variational basis is. The efficiency of the 
variational procedure is manifest in the overlap of the vacuum with the lowest-lying stringy state (i.e. the equivalent 
of the coefficient |coP in Eq. (HH])), which would be one in the ideal case in which the variational procedure identified 
the exact state we are interested in. In the considered fuzzying scheme, with a careful choice of the parameters a and 
6 it is possible to reach overlaps of the order of 0.9-0.95, which makes the contribution of excited states negligible 
already at distances of the order of two lattice spacings. Physically, the process of blocking and smearing allows us 
to build variational trial states on the scale of physical distances, while simple Polyakov loop correlators probe the 
physics on the scale of the lattice spacing, which is sensitive to ultraviolet fluctuations. 
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FIG. 5: Spatial and temporal string tensions at various lattice sizes as a function of amo. Also shown is the infinite volume 
estimate. 



We remind that we use Eq. ([5^ for extracting the string tension from Eq ^. Our results for the string tension 
are plotted in Fig. [SJ For the smearing and blocking procedures, we have used a — 0.4 and S = 0.16. As expected, 
at high mass the spatial and the temporal string tensions agree and are independent of the lattice size. The signal 
provided by the Polyakov loop correlator is clean and the overlap |cop is of order 0.9 — 0.95. As amo is reduced, on 
the smaller lattices the spatial and the temporal string tensions depart, and the overlap becomes of order 0.8. The 
disagreement between the temporal and the spatial string tensions is an indication that finite size effects are starting 
to play a major role. As a matter of fact, if we simulate at fixed values of aniQ on increasingly larger lattices, we find 
that the spatial and the temporal string tensions eventually agree, and that they also agree with the temporal string 
tension of the smaller lattices. This shows that the temporal string tension is less affected by finite size effects. This 
is hardly surprising, since the onset of the departure between the spatial and temporal string tension arises where the 
system goes from the S- to the A-phase. This also confirms that in the A-regime strictly speaking it is not correct 
to talk about a spatial string tension, since the mass of the Polyakov loop is not associated to confining strings. 
However, for convenience we shall still use the string language. For a given volume, if the mass is lowered below the 
onset of spatial deconfinement, the agreement of the temporal string tension and the string tensions computed on 
larger volumes is lost, and the former flattens out. This kind of finite size effects appears when the correlation length 
associated to the string tension becomes of the order of the spatial lattice size. In fact, for the plateau values we find 
a^at ~ 4:N^'^ = ^s^^- Taken at face value, this would imply that in the thermodynamic and chiral limit Uy/ot — 0. 
However, since these results have been obtained in a phase were finite size artefacts play a major role, a confirmation 
of this statement on larger lattices is necessary before we can conclude that there is no asymptotic string tension in 
the massless limit. 

We estimate the string tension at infinite volume by choosing the determinations that are reasonably safe from finite 
volume effects. Our strategy is based on the following observations. At fixed mass and in the infinite volume limit, 
the temporal and spatial string tensions must coincide. At fixed volume we observe that at large enough mass, the 
temporal and spatial string tensions coincide. In general the spatial string tension can be determined more accurately, 
since correlators of shorter Polyakov loops have smaller relative errors. On the other hand, the temporal string tension 
is less affected by finite volume effects and in particular is always well defined. Whenever the spatial and temporal 



^ Alternatively, Eq. JSTJ can be used, since it gives results that are compatible well within errors with Eq. |(52j; on the contrary, Eq. II50I I 
gives discrepancies of up to 20% for loops winding the spatial directions. 
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string tensions agree within one standard deviation, we consider the weighted average of the two: 

I \/57 

^^(A^F_^IW^ (56) 

AV^ = ^ ^ ^ = . (57) 

In the A-phase, where the temporal and spatial string tensions do not agree anymore, the correct string tension to be 
considered is the temporal one. At those values of the mass, for which more than one volume is available, the result 
on the largest volume has been considered. 

More in detail, at the mass amo = —1 both the 24 x 12^ and 16 x 8"^ lattices are available, and the two temporal 
string tensions (plus the spatial string tension on the larger volume) are in agreement at the la level. For this mass 
and larger ones we therefore expect that the temporal string tension, as determined on the 24 x 12"^ and 16 x 8"^ 
lattices, is affected by smaller finite-volume effects than the statistical indetermination. In particular, the string 
tensions on the 24 x 12'^ lattice at masses amo = ^1 a-nd —0.95, and on the 16 x 8^ lattice at masses amo = 
—0.975, —0.9, —0.75, —0.5, —0.25, 0, 0.25 are good estimates for the infinite volume limit. 

Moreover, at the mass amo = —1.15 both the 32 x 16'^ and 24 x 12^^ lattices are available, and the two temporal 
string tensions agree at the la level. Repeating the argument above, the string tensions on the 32 x 16'' lattice at 
mass amo = —1.15, and on the 24 x 12'^ lattice at masses amo = —1.125, —1.1, —1.075, —1.05 are good estimates of the 
infinite volume limit. Below mass amo = —1.15 we have no control on the finite volume effects, so we just discarded 
those points. 

The choices described above, and the results for the estimate of the string tension at infinite volume are summarized 
in Table ITVl All the measured string tensions and the infinite volume estimates are also plotted in Fig. [5] 
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w.a. on 24 x 12^ (BO) data 


24 X 12^ 


0.975 


0.2066(97) 


temporal s.t. on 16 x 8^ (Al) 




1 


0.1851(33) 


w.a. on 24 x 12^ (Bl) data 


24 X 12^ 


1.05 


0.1587(23) 


w.a. on 24 x 12^ (B2) data 


24 X 12^ 


1.075 


0.1455(19) 


w.a. on 24 x 12^ (B3) data 


24 X 12^ 


1.1 


0.1205(56) 


temporal s.t. on 24 x 12^ (B4) 


24 X 12^ 


1.125 


0.1130(54) 


temporal s.t. on 24 x 12^ (B5) 


24 X 24 X 12^ 


1.15 


0.0790(34) 


temporal s.t. on 32 x 16^ (CO) 


32 X 16^ 



TABLE IV: This table shows the infinite volume estimate for the string tension (2nd column). Where different volumes are 
available for the same bare mass, the larger one has been chosen. If the temporal and spatial string tensions agree within 
one standard deviation, a weighted average (w.a.) between the two has been chosen, as explained in the text. Otherwise the 
temporal string tension (s.t.) has been chosen. For the bare mass amo ~ —1, the temporal string tensions on the 16 x 8'^ and 
24 X 12'' lattices agree at the la level: this suggests that for all the string tensions computed at this and higher masses on 
the 16 X 8^ lattice, finite volume effects are negligible. Also the temporal string tension computed on the 32 x 16'' lattice at 
amo = —1.15 coincides at la with the one computed on the 24 x 12^ lattice. 



VI. STATIC FORCE AND POTENTIAL FROM WILSON LOOPS 

A different way to compute the string tension is via the expectation values of Wilson loops. The advantage of this 
method is that it can show whether the mass extracted via Polyakov loop correlators is in fact related to the existence 
of an asymptotic string tension at fixed lattice geometry. As we will see, the disadvantage is that the numerical results 
for the string tension extracted from Wilson loops generally have larger statistical errors. 
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We consider the generic off-axis Wilson loop W{T,'R), describing a quark-antiquark pair separated by a spatial 
distance R = {X, Y, Z) and propagating in straight line in the temporal direction. It will be useful to schematically 
decompose the close parallel transport in its spatial Us[t;r ^ r + R] and temporal Ut[t ^ t + T; r] components: 

W{T,R) = TT{Ut[t t + T;r]''Us[t]r ^ r + Ii]Ut[t ^ < + T; r + R]C/^[f + T; r ^ r + R]^} . (58) 

The off-axis component Us[t;r — >■ r -|- R] is computed following the algorithm in Ref. [80|. In order to reduce the 
noise in the static potential, we build the Wilson loops with smeared link variables. We choose a single step of HYP 
smearing [8l|; the smeared link variable is a function of all the links belonging to the unit hypercubes adjacent to 
the original link. We found the HYP smearing effective enough for our purposes, even not implementing a variational 
method. The main disadvantage is that it deforms the static potential at short distances; however it does not affect 
the determination of the string tension, which is our main goal. 

In the Hamiltonian gauge ((7(x, 0) = 1), the expectation value of WiT^TL) becomes the time-correlator for the 
operator M(R) — Q{r)Us[r — ?> r + R](5(r + R) which creates a heavy quark-antiquark pair connected by a string: 

(T^(T,R)) = ^^Tr[e-(-^'-^)^Af(R)te-^^A/(R)] = 
Z{Lt) 

- ^Y.\(^,R\Mmn^}\'e-^^^-e-^^''-^^^^-^-^ , (59) 

where we have inserted the gauge invariant states |m) with energy i?,„, and the states in presence of a quark-antiquark 
pair |n, R) with energy Vn{R)- The energy Vo(-R) of the fundamental state in presence of the quark-antiquark pair is 
the static potential. Z{Lt) is the partition function, which can be decomposed in terms of the gauge invariant states 
\m): 

Z{Lt) = Tre-^'" = e"^'-^-" . (60) 

m 

If the temporal extension Lt is large enough, we can identify three regimes for the Wilson loops while T changes: 

1. Small values of T with respect to the temporal extension (T <ti Lt). This is the usual zero temperature limit. 
In this case only the vacuum survives among the |m) states: 

{WiT, R)) ^ J2 |an(^)l'e-^[^"(«)-^°l , with ^ |a„(i?)p - 2 , (61) 

n n 

where the relationship for the coefficients Q!„ comes from (14^(0, R)) = 2. It is interesting to notice that, since 
{W{T, R)) < 2 for each value of T, Vn{R) must be larger that the vacuum energy Eq for each value of i? > 
(because of the lattice discretization, the potential is bounded from below). Therefore in this regime, the Wilson 
loop is decreasing in T at every fixed R. 

2. Values of T comparable with the temporal extension {Lt — T <C Lt). In this case only the |0, R) survives among 
the states with external charges: 

(I^(T,R)) ~ e-^*[^o(«)-^o] ^ |^^(^)|2g-(L,-r)[s„-v„(fl)] (g2) 

m 

Assume that the sum is dominated by some to. If the quantity Em — Vq{R) is positive, the Wilson loop at fixed 
R is increasing as T approaches Lt, towards its extremal value (VF(L(, R)) ~ J2m \Pm{Li)\'^ < 2. If the quantity 
Em — Vo{R) is negative, the Wilson loop at fixed R is decreasing towards its extremal value. Notice that, since 
the states propagating around the torus arc not the same as the ones propagating inside the Wilson loop, the 
Wilson loop is not symmetric in T around Lt/2 as it usually happens for other correlators. 

3. Intermediate values of T (a ^ T <C Lt). In this case the Wilson loop reduces to a single exponential: 

(I^(T,R)) ~ A(R)e-'^[^''(-«)-^''l . (63) 

This is the useful regime which we will try to identify in our numerical simulations to extract the static potential. 

The first and second regions are always visible if the temporal length T of the Wilson loop is too small or too large. 
With the extra difficulty that the Wilson loops are not symmetric in T — >■ Lt — T, the computational problem for 
the static potential is similar to the one for other correlators: it is important to have a large enough lattice in such a 
way that the third region opens up in the middle. An effective potential is therefore defined, with the property that 
it shows a plateau in the third region (if visible), and the value of the plateau is actually the static potential. We use 
two different methods to extract the effective potential from the Wilson loops. 
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Potentiall . The easiest method consists in defining an effective potential as: 

If we can see a plateau in the effective potential as a function of T, it means that we can isolate the single- 
exponential region. The value V^(R) is then extracted by fitting the plateau of the effective potential with a 
constant. Notice that since the lattice breaks rotational invariance, we consider the potential as a function of 
R, and not of its module only. An unbiased estimate for the average of the potential, and an estimate of its 
error are obtained by applying Eq. to a set of bootstrap ensembles. 

Potential2 . When the single-exponential region is not visible, we use the Prony's method for taking into account 
also the first excited state. We refer to the literature for the general idea, while we summarize here the used 
formulae. Having chosen a value of R, for every value of T we solve the following second-order equation: 

[W{T, R)W{T + 2a, R) - W{T + a, R)^] (65) 
+ [W{T + a, Ii)W{T + 2a, R) - W{T, R)W(T + 3a, R)] x+ (66) 
+ [W{T + a, R)VK(T + 3a, R) - W{T + 2a, Rf] = . (67) 

If xq is the largest solution (but smaller than 1), we define an effective potential as 

Feff(T,R) --ilogcco . (68) 
a 

As for the previous method, the value 1^(R) is extracted by fitting the plateau of the effective potential with 
a constant. The whole procedure is implemented via a bootstrap, in order to get an unbiased estimate for the 
average of the energies, and an estimate of its error. 

Whenever we can compute the potential with both methods we observe that they always give compatible results, 
but slightly smaller errors and better determinations of the plateaux are obtained with the method Potential2. Even 
though the string tension can be in principle extracted from a linear fit of the static potential in a large distance 
region, it is instructive to determine it also from the force F = — At large R, F ^ a. For the determination of F, 
we use the following methods. 

Forcel . We use generalized Creutz ratios to define an effective force with off-axis Wilson loops: 

rrr. N 1 , (W {T + a, R + an)) (W {T ,R)) 

We use n vectors of the form (1,0,0), (1, 1,0) and (1, 1, 1) and permutations. We identify the plateau of the 
effective force as a function of T, and we fit it with a constant F{Rj). The statistical error is determined by a 
bootstrap procedure. The improved distance Rj is defined as in Ref. [s^ to be: 



Ri 



^^ G(R + n)-G(R) 



-1/2 

(70) 



where G(R) is the three-dimensional free-scalar propagator on the lattice. 

Force2 . Plateaux in the Creutz ratios are visible only in a region where only a single exponent dominates in the 
expansion of the Wilson loop. In most of the cases we need to take into account the first excited state. An 
effective force can be defined by using the effective potentials computed by the method Potential2: 

p (T-n Kff(T,R + an)-Kff(T,R) 

Feff(r,R,an) = . (71) 

a|n| 

We identify the plateau of the effective force as a function of T, and we fit it with a constant F{Ri). Expectation 
value and error of F{Ri) are estimated by means of a bootstrap procedure. 

The static potentials presented in this section have been computed with the method Potential2. In Fig. |6l two 
typical effective potentials are shown, together with the fit range and the result of the constant fit. 

The static potentials for all the simulations listed in the last column of Table IIVI are plotted in Fig. [71 The 
corresponding forces are separately plotted in Fig. [H Although in principle the static potential or the force can 
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FIG. 6: Effective potentials at R = (4a, 0, 0) and R = (6a, 6a, 6a) computed on the 32 x 16'' lattice at bare mass amo = —1.15. 
Plateaux have been chosen in the white regions. The red lines represent the values of the potential (with errors), obtained by 
fitting the effective potential with a constant in the white region. Errors have been computed with a bootstrap procedure. 



be used to extract the string tension, in practice a reliable result can not be obtained from those quantities, the 
most likely explanation being either that our data are not accurate enough or that the lattice sizes explored are too 
small for the plateau to be free from systematic errors. A variational procedure like the one used for extracting the 
string tension from correlators of Polyakov loops (described in Sec. |V| might be helpful also for the static potential 
computation. Although we are unable to perform a comparison between the string tensions extracted with the two 
different methods, we can still check that the static potentials and the forces are compatible with the string tensions 
reported in Table HVl 




FIG. 7; Static potentials computed with the method Potential2 (see Sec. IVIfl from Wilson loops with HYP smearing. Data 
for bare masses —0.95, —1, —1.05, —1.075, —1.1 are obtained on a 24 x 12'^ lattice; data for bare mass —1.125 are obtained 
on a 24 X 24 X 12^ lattice; data for bare mass —1.15 are obtained on a 32 x 16^ lattice. Superimposed is the function 
V{R) = crR + fi + c/R, where a has been taken from Table HVl (the two curves correspond to cr — Act and a + Act), while fi and 
c have been obtained with a fit in the region R > 3a. 
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We fit the static potential with the function: 

V{R) =<jR + fi+-^ , (72) 
-ft 

assuming the string tensions shown in Table IIVI in the range R > 3a. The results of the fits are shown in Fig. [T] 
Since we do not want to assume at this stage a particular effective string theory (and anyhow the HYP smearing 
introduces spurious 1/R effects) , the coefficient of the l/R term becomes an extra parameter in the fitting procedure. 
In all the cases we have investigated, the string tension computed via Polyakov loop correlators captures correctly the 
large distance behavior of the static potential. 

In Fig. m the forces are plotted together with the values of the string tension. Although the errors on the force 
are in some cases quite large and only qualitative statements are possible, it can be seen also in this case that the 
string tension computed via Polyakov loop correlators always captures the large distance behavior of the force itself. 
The force always shows a plateaux at large distances, with a central value often in striking agreement with the string 
tension computed from Polyakov loop correlators. This might indicate that our analysis overestimates the statistical 
errors. 

VII. GLUEBALLS 

Glueball masses are extracted from a variational procedure similar to the one used for Polyakov loops and based on 
the same fuzzying scheme. At the link level, we consider a collection of closed elementary loops transforming according 
to the irreducible representations of the symmetry group of the cube, to which rotational symmetry is broken on the 
lattice (s^ . The variational procedure is then built by replacing the original links with those obtained after smearing 
and blocking. Then, in each channel a matrix of connected correlators is constructed, whose eigenvectors with the 
highest eigenvalues are almost the pure eigenstates of the Hamiltonian with the lowest masses. We have implemented 
this technique using as starting operators the plaquette and the length-six planar closed contour to build the A, E and 
T irreducible representations of the cubic group. The lowest-lying state in the A channel corresponds to the lightest 
0"'""'" glueball in the continuum limit, while both the E and the T lightest states give the lightest 2++ glueball mass 
in the continuum limit. 

Our results for the 0++ glueball are illustrated in Fig. [HI We discuss only the results on the 16 x 8'^ and the 24 x 12'^ 
lattices, since on the 32 x 16^ lattice for all values of the simulated bare fermion masses our system is in the A-phase. 
As it can be seen from the plot, as amo is decreased, in the S-phase aMo++ monotonically decreases. When the 
system starts to develop double peaks for the Polyakov loop distribution in a spatial direction, the mass of the 0++ 
glueball first goes to a short plateau and then moderately increases, to drop dramatically to much lower values at 
lower amo. The overlap between the masses measured on the two lattices at otoq = —0.95 suggests that as long as we 
stay in the symmetric phase, finite size effects are under control. Hence, as our best estimate for the infinite volume 
limit of aMo++, we take the values on the 16 x 8'^ lattice for amo > —0.95 and the values on the 24 x 12'^ lattice 
for —1.05 < amo ^ —0.95. Since for amo < —1.05 in both cases our system is in the A-phase, we do not take into 
account the corresponding values of amo+-f in the following analysis. 

The mass of the 2++ glueball as a function of amo is shown in Fig. [TOl These results have been obtained using 
operators transforming according to the E representation of the cubic group. Our analysis in the T channel gives 
compatible results in all cases. Our data show that the 2++ glueball is heavier than the 0"'"+ in the symmetric phase, 
but dramatically decreases to very low values of the mass (well below the mass of the 0++) at the onset of the A-phase. 
Deeper in the A-phase, the two states appear to be degenerate. As the figure shows, no good control over finite size 
effects can be reached on our lattices for the 2++ mass. For the sake of completeness, we still provide an estimate for 
its mass at infinite volume, but this is likely to be quite rough. Hence, the 2++ glueball will play a marginal role in 
the interpretation of our results. 

Our numerical estimates of aMo++ and aM2++ in the infinite volume limit are reported in Tab. [V] The degeneracy 
between the two states at m = —1.05 together with the impossibility of establishing whether the system is in the 
S-phase (see Tab. [TTI would suggest to disregard glueball masses at this value of the bare mass. However, since 
this point was part of our analysis in Ref. [5a], where the lower statistics masked the issue, in order to facilitate a 
comparison with our previous work, we chose to keep it also in our current analysis. The reader should bear this in 
mind for the discussion of our results. 

VIII. HYPERSCALING SCENARIO AND LOCKING SCALE 

Our infinite volume estimates for the PS (at each value of the bare mass, we choose the PS mass computed on 
the largest volume in [131), the O"*""*" and the 2++ glueball masses and a^/"^ as a function of the PC AC mass am (see 
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—amo 


aMf,++ 


aM2++ 


- 0.25 












0.5 


1.045(70) 


1.93(19) 


0.75 


0.919(63) 


1.27(12) 


0.9 


0.666(44) 


0.874(77) 


0.95 


0.793(41) 


1.129(85) 


1 


0.658(56) 


0.886(73) 


1.05 


0.510(33) 


0.513(35) 



TABLE V: Infinite volume estimates of aMQ++ and 0M2++ . Values extracted on a 16 x 8'' lattice have been used for amo > —0.95 
and values extracted on a 24 x 12"^ lattice for amo < —0.95. 

Refs. [1^ [13 for a definition of this quantity) are reported in Fig. [TTJ As noticed in Ref. [56^, the data show a clear 
hierarchy in the spectrum, with the mesonic scale well above the gluonic scale. Since over the range of investigated 
masses acr^/^ changes by a factor of five, the effect of the fermion determinant as the mass is decreased is an essential 
component of the dynamics in this theory. Hence, the simple quenched scenario, according to which the theory would 
be QCD-like and the hierarchy in the spectrum is due to large fermion masses, can be excluded. In fact, the spectrum 
looks similar to the hyperscaling scenario at high locking mass Afiock sketched in Fig. [T] (right). In this section we 
shall show that indeed that scenario provides the right description of the spectrum of this theory. 

Up to subleading corrections, the hyperscaling scenario implies the independence of ratios of physical quantities 
from the fermion mass in the scaling region. In Fig. [12] we plot the ratio Mps/cr^^^- This quantity shows a plateau 
Mps/a^^"^ ~ 7.5 for aMps < 1.25, supporting the idea that gluonic and fermionic masses are not parametrically 
independent in this region but are both proportional to the RG-invariant fermion mass M (again, this is what we are 
calling locking) . The independence from M of the ratio of spectral quantities in the scaling region can be observed 
also in the ratio My/Mpg (Fig. [T5)) . Again a plateau develops for aMpg < 1.25. The existence of these plateaux is a 
clear indication of the spectrum behaving as predicted by hyperscaling. Moreover, the value of both Mps/a^^^ and 
My/Mps in the scaling region suggests that the locking scale Miock is large. If this is the case, SU(2) gauge theory 
with two adjoint Dirac fermions should look like an heavy fermion system for all values of the fermion mass. 

In order to verify this scenario, we can compare our dynamical results with results obtained in the quenched theory. 
In this process, a crucial point is to match properly the bare parameters in the two theories, since the physics of the 
dynamical and quenched systems is expected to be equivalent once the physical scale is matched (see Eq. (15^ ). For 
this reason, we need to tune the bare parameters of the quenched simulations, namely the gauge coupling /S^'^^ , and 
the mass of the valence fermion am'^\ so that we can match two independent quantities between the two theories. 
For this matching, it is convenient to choose two observables with a sharp dependency on each of the bare parameters, 
in order to achieve the best possible tuning between the two theories. In our study, we have required the quantities 
aa^^"^ and aAIps to be equal in the dynamical and in the quenched systems. Given that in the dynamical theory with 
adjoint fermions the fundamental string cannot break, this quantity is a natural candidate to fix the coupling in 

the quenched simulation. As far as the valence quark mass am'^^ is concerned, we chose to determine the quenched 
value by fixing the pseudoscalar mass, because in the fermionic sector this is the quantity we have the best control 
upon. 

The procedure to compare the quenched and dynamical theories requires the following steps: 

1. find the value Z?*^'-' of the coupling for the quenched theory, in such a way that the string tension in lattice units 
matches the dynamical value; 

2. find the value arn^^ of the bare fermion mass for the quenched theory, in such a way that the PS meson mass 
matches the dynamical value aMps; 

3. compute the rest of the spectrum, for example the ratio of the PS and V meson masses or the glueball masses, 
in the dynamical theory with bare parameters (/3, amo) and in the quenched theory with bare parameters 
(/3^*\ ami,''''), and compare the results. 

In practice the program as outlined above requires a fine tuning of the bare parameters, and turns out to be a 
highly expensive computational task. Instead of an exact matching of the parameters, we have performed a scan in 
the parameters of the quenched theory. The lattice parameters at which quenched simulations have been performed 
are reported in Tab. IVII The range of the scan of /3'^*-' is chosen in order to include all the string tension values of our 
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FIG. 8: Forces computed with the method Force2 (see Sec. IVIf) . Superimposed is the string tension from Table [TVl 
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FIG. 9: The mass of the 0"*"^ glueball in lattice units, 
aAlQ++ , measured at various values of bare quark mass 
amo on a 16 X 8^ and on a 24 x 12^ lattice. 



FIG. 10: The mass of the 2^'^ glueball in lattice units, 
aM2++, measured at various values of bare quark mass 
amo on a 16 X 8^ and on a 24 x 12^ lattice. 
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FIG. 11: The spectrum of the theory as a function of the PCAC mass am. The mass of the vector is not shown, since on the 
scale of the figure this state appears to be degenerate with the PS. 



dynamical simulations. The upper bound of this window is simply given by /3 = 2.25, the value of /3 for our dynamical 
simulations. This is a consequence of the string tension being an increasing function of the bare fermion mass, so an 
infinite mass simulation corresponds to a string tension of a pure gauge system at the same /3 = 2.25. 

For each choice of /^^'^ we measure the string tension and the 0++ and 2++ glueball masses (following e.g. Ref. 0]). 
Quenched glueball masses have been interpolated using the ansatz 

M^=Ao + A,a^a, (73) 

with Aq and Ai respectively the leading (constant) and subleading (©(a^)) coefficients in the extrapolation to the 
continuum limit. 

On the same gauge configurations, we measure the quenched PS mass and the My/Mps ratio, for a set of values of 
arrS^^ covering the entire interval of PS masses appearing in the dynamical calculation. We then create an interpolating 
function for the central value of the ratio My/Mpg. 

To obtain an error on this estimate, we create two other interpolating functions for the maximal and the minimal 
value of the quenched estimate My /Mps set by the statistical error, so that for each choice of the pair {Mps, a), we 
can read the corresponding range of values for My/Mps. To take into account the indetermination in our estimate 
of aMpg and a^a, we consider a region within one sigma around the central value for those quantities: in this region 
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FIG. 12: The ratio Mps/(J ' as a function of Afr 



the difference between the maximum value of the maximal interpolating function and the minimum of the minimal 
interpolating function provides us with an estimate for the error on My/Mps in the quenched theory. 

By means of the interpolating functions, we can read the value of Alc/y/a and My/Mps in the quenched simula- 
tions at values of ay/a and aMpg (the latter being relevant only for the My /Mps ratio) obtained in the dynamical 
simulations. 

In Tab. IVIII we report the values of the mass of the glueballs for the quenched theory at the values of the string 
tension obtained in the dynamical theory. In Fig. [TJ] we show the comparison between the dynamical glueball values 
and the interpolating functions obtained from the quenched theory. Except for the last point, for which, as discussed 
in Sect. IVIII the dynamical simulations are probably in the A-phase, the agreement between the quenched and 
dynamical spectra at the same physical scale (in units of the ultraviolet cutoff) is striking. This supports the idea 
that the low-energy dynamics of the theory with dynamical fermions is well described by a pure Yang-Mills theory, 
which is evidence for a locking mechanism with a large Afiock taking place in this theory. 

The relevant interpolated quenched results for the ratio My/Mps as a function of Afps are reported in Tab. IVIlTl 
and compared with the dynamical results in Fig. 1131 The quenched and dynamical data have a remarkable overlap for 
all the points except for the last one, where strong finite-size effects are expected to affect both the quenched and the 
dynamical simulations. Together with the plateau in the ratio developing for aAfps < 1.25, the comparison confirms 
once again a locking mechanism taking place at large Miock- 



IX. THE CHIRAL CONDENSATE ANOMALOUS DIMENSION 



The hyperscaling scenario supported by our data seems to imply the existence of an infrared fixed point. However, 
since the evidence for the hyperscaling and the locking of the mesonic and gluonic spectra is still over a small range 
of am, simulations at smaller masses and large volumes are needed to confirm the trend identified so far. 

If the theory is IR conformal, all the spectral quantities scale as for a unique value ofp=l/(l-f7:„), with 7* 
the anomalous dimension of the condensate. Hence, in this case 7* is physically well-defined. 

In order to build phenomenologically viable Technicolor models, a large anomalous dimension is generally required. 
From a purely theoretical point of view < 7* < 2, where 7* = corresponds to the non- interacting case and 7* = 2 is 
the bound imposed by unitarity; a value 7^ « 1 might reconcile Technicolor with high-precision data for the Standard 
Model. The determination of 7* is then one of the goals of lattice simulations of BSM strong dynamics. 

In order for us to be able to extract a scaling exponent, the simulations must be performed in a region of sufficiently 
small masses. The exact extent of the scaling region (which also depends on the observable being analyzed) is only 
known a posteriori. On the lattice, the problem is complicated by the explicit breaking of conformal invariance due 
to the finite size of the system. This can however be turned into a powerful tool for determining the exponent of the 
scaling with the mass using a technique commonly known in Statistical Mechanics as Finite Size Scaling (FSS). FSS 
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FIG. 13: Comparison of the ratio AIv/Mps as function of aMps in the quenched and the dynamical theory. The dynamical 
ratio approaches the infinite-mass value of one at large aAIps and develops a plateau (signaling hyperscaling) at aAfps < 1.25. 
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FIG. 14: Comparison of the masses of glueballs of the dynamical theory with the interpolating functions of the quenched theory. 
Among the plotted data, the dynamical points at aa^''^ = 0.2405(58) and aa^^^ = 0.1587(23) correspond to the lightest masses 
respectively on the 16 x 8'^ and 24 x 12^ lattices. Hence, the glueball masses for those string tensions are affected by the largest 
finite-size effects. 



states that the dimension of the system is a relevant scaling variable with mass dimension -1. Hence, the asymptotic 
scaling formula 

aMxoc{amy, p=l/{l + j,), (74) 

where Mx is a spectral quantity of the system, on a finite lattice of spatial extension Ls = aNg and in the regime 
is — >■ oo and m — becomes 

MxLs = f{x), x = Ns{am)P , (75) 

i.e. the product MxLg is a universal function of the scaling variable x. A simple consequence is that the ratio of two 
spectral quantities is expected to be a universal function of MxLs for any spectral quantity Mx- Note that this is 
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V 




2.25 


0.4231(25) 


24 X 12^ 


-1 65 -1 6 -1 55 -1 5 -1 45 -1 4 -1 35 -1 3 
-1.25,-1.2,-1.175,-1.15,-1.125,-1.1,-1-075, 
-1 05 -1 025 -1 -0 75 -0 5 -0 25 25 5 


2.4265 


0.2388(9) 


24 X 12^ 


-1.2,-1.175,-1.15,-1.125,-1.1,-1.075,-1.05, 
-1 025 -1 -0 75 -0 5 -0 25 25 5 


2.5115 


0.1768(8) 


24 X 12^ 


-1.2,-1.175,-1.15,-1.125,-1.1,-1.075,-1.05, 
-1.025,-1.0,-0.75,-0.5,-0.25,0.0,0.25,0.5 


2.6 


0.13395(62) 24 x 12^ 


-1.2,-1.175,-1.15,-1.125,-1.1,-1.075,-1.05, 
-1.025,-1.0,-0.75,-0.5,-0.25,0.0,0.25,0.5 


2.62 


0.1258(7) 


24 X 12=* 


-1.3,-1.25,-1.2,-1.15,-1.1,-1.05,-1.0,-0.75, 
-0.5,-0.25,0.0,0.25,0.5 


2.68 


0.1035(7) 


32 X 32=* 


-1.25,-1.2,-1.15,-1.1,-1.05,-1,-0.95,-0.9, 
-0.85,-0.8,-0.75,-0.5 



TABLE VI: Bare parameters and volumes used for quenched simulations. For each /S*^'-' we report also the measured string 
tension. 





0+ + 




0.4053(51) 


1.30(4) 


2.13(15) 


0.352(11) 


1.18(5) 


1.89(14) 


0.3561(42) 


1.185(35) 


1.90(11) 


0.2794(63) 


0.975(35) 


1.53(8) 


0.2405(58) 


0.855(35) 


1.325(65) 


0.2184(23) 


0.785(15) 


1.21(4) 


0.1851(33) 


0.675(25) 


1.025(35) 


0.1587(23) 


0.585(15) 


0.885(25) 



TABLE VII: The results of the interpolation procedure for the quenched data of the mass of the glueballs at the values of aa^' 
equal to the ones found in the dynamical simulations. The errors on the third column have been obtained from the variation 
of the numerical results on the maximal and the minimal interpolating function. 
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true for both the S and the A-phases of the system. In Fig. [12] we show the ratio ct/ /Mps as a function of MpsLg. 
The universahty of the ratio is verified up to values of MpsLs — 12. 

If the system is in the scaling region (for which we have support from our data) and the infinite-volume estimates 
for spectral quantities are correct, Eq. (|74p can be used for determining 7,. We have shown in Sect. |V] that our 
determination of the string tension is reasonably under control. Hence, we perform a fit of our data for cr^/^ using 
Eq. (fM]). With a good quality of the fit (see Fig. [TH]), we find 7* = 0.22(6). Both horizontal and vertical data errors 
have been taken into account by implementing a bootstrap procedure. The fit has been performed on the lightest four 
points, and then progressively increasing the fitting region to include ten points (all but the last three in Fig. 1161) . 
No systematic trend has been observed when enlarging the fitting region. The quoted value of 7* is a conservative 
estimate compatible with all the values obtained using the fitting procedure described above. 

The value we find for 7^, is compatible with determinations obtained in the same theory using other quantities (e.g. 
related to mesonic physics [53) or independent techniques like the Schroedinger functional t28|]. These results clearly 
favor the existence of a genuine IR fixed point for this theory. The fact that independent measurements of 7, fall all 
in the same window of values is a clear message for model building. 

X. CONCLUSIONS 

In this work, using numerical simulations of the lattice model for several sizes of the system and a wide range 
of fermion masses, we have shown that at sufficiently low masses the spectrum of Minimal Walking Technicolor is 
consistent with the existence of an infrared fixed point. In particular, for this specific realization of locking, the theory 
at large distances is isospectral to a Yang-Mills SU(2) theory, where the dynamically generated scale of the pure gauge 
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0.2066(97) 


1.4094(42) 


1.0389(40) 


0.1851(33) 


1.3493(28) 


1.0379(21) 


0.1587(23) 


1.1874(28) 


1.0419(37) 


0.1455(19) 


1.0811(31) 


1.0464(30) 


0.1205(56) 


0.9613(35) 


1.0396(70) 


0.1130(54) 


0.8017(41) 


1.053(13) 



TABLE VIII: The results of the interpolation procedure for the quenched data of M^^ /Mp'g at the values of aa^^^ and aMps 
equal to the ones found in the dynamical simulations. The errors on the third column have been obtained from the variation 
of the numerical results on the maximal and the minimal interpolating function. 




FIG. 15: The ratio ay^/Mps as a function of Mpsl/a at FIG. 16: aa^^^ as a function of am. A fit of the data 

various lattice sizes. In the shaded region the crossover to Eq. (|74p is also shown. In particular, the two lines 

between the S and the A phase takes place. represent the extremal values 7* — 0.16 and 7* — 0.28. 



theory is determined by the fermion mass in MWT, and turns out to be smaller than the latter. To confirm this 
scenario would require to extend our study to much smaller fermion masses, down to values that are not accessible 
at present to lattice simulations. Another technical limitation of our study is the simulation at a fixed value of /?: in 
order to verify that lattice artefacts are not distorting the physical picture, further studies closer to the continuum 
limit should be performed. 

Assuming the existence of the IR fixed point, we have determined the anomalous dimension of the condensate, which 
is found to be 7* = 0.22(6). This value is in agreement with other independent determinations, which strengthen the 
conclusions that the theory is infrared conformal. The value of 7* for this theory is probably too small for conventional 
Technicolor scenarios, although alternative scenarios compatible with a small anomalous dimension can be devised 
(see e.g. [1^1 )• It would be desirable to include larger lattices in our FSS analysis. 

Finally, we notice that a FSS analysis performed in a SU(3) gauge theory with two fermion flavors in the two-index 
symmetric representation (which for SU(2) coincides with the adjoint representation) finds 7* ~ 0.5 [30] . Assuming 
that, as stated in Refs. [Ij, [IJ, that this theory is infrared conformal (however, see Refs. (13, HH] for alternative 
scenarios), this might imply that 7* for two- index symmetric fermions is an increasing function of the number of 
colors N. If this is the case, it would be interesting to determine whether 7* becomes of order one for large enough 
values of N and whether 7* also increases with N for adjoint fermions. 
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